Circuit theory for decoherence in superconducting charge qubits 



IBM T. J. 



Guido Burkard 

Watson Research Center, P. O. Box 218, Yorktown Heights, NY 10598, USA 



in ■ 

o : 
o . 

(N ■ 

: 

<: 

oo ■ 

(N : 



X3 

o 

> 
OO 

oo 

00 

o 

O 



I 

O 

o 



X 



Based on a network graph analysis of the underlying circuit, a quantum theory of arbitrary 
superconducting charge qubits is derived. Describing the dissipative elements of the circuit with 
a Caldeira-Leggett model, we calculate the decoherence and leakage rates of a charge qubit. The 
analysis includes decoherence due to a dissipative circuit element such as a voltage source or the 
quasiparticle resistances of the Josephson junctions in the circuit. The theory presented here is dual 
to the quantum circuit theory for superconducting flux qubits. In contrast to spin-boson models, 
the full Hilbert space structure of the qubit and its coupling to the dissipative environment is taken 
into account. Moreover, both self and mutual inductances of the circuit are fully included. 



I. INTRODUCTION 



Various types of quantum bits with Josephson june- 
tions in superconducting (SC) circuits are now inves- 
tigated in theoretical and experimental studies P, 
The two types of macroscopic SC qubits, the charge 
^, "S", "^,0 13 and flux d lis m m qubits, are dis- 
tinguished by the relative size of the charging energy Ec 
and the Josephson energy Ej of their junctions [13 ■ In 
flux qubits, also known as persistent-current qubits, the 
Josephson energy dominates, Ej ^ Eq, and the state 
of the qubit is represented as the orientation of a per- 
sistent current in a SC loop 0, 0, 0|- In con- 
trast to flux qubits, charge qubits operate in the regime 
Ec ^ Ej, and are represented as the charge state of a 
small SC island (presence, or absence, |0), of an extra 
Cooper pair) which is capacitively coupled to SC leads 
|3','?,"5",'a,|l|| (Fig.d. The quantronium is a charge 
qubit that operates in a regime close to Ec ~ Ej. 

Both types of SC qubits suffer from decoherence that is 
caused by a several sources. In flux qubits, the Johnson- 
Nyquist noise from lossy circuit elements (e.g., current 
sources) has been identified as one important cause of 
decoherence 0, 0, . A systematic theory of deco- 
herence of a qubit from such dissipative elements, based 
on the network graph analysis [21| of the underlying SC 
circuit, was developed for SC flux qubits 0, and suc- 
cessfully applied to study the effect of asymmetries in 
a persistent-current qubit p^ . Decoherence in charge 
qubits has previously been investigated using the spin- 
boson model P, 113 ■ 

Here, we develop a general network graph theory for 
charge qubits and give examples for its application. As 
in the case of the circuit theory for flux qubits, we are 
not restricted to a Hilbert space of the SC device which is 
a priori truncated to two levels only. In other words, in 
contrast to the spin-boson model, our theory is capable of 
describing leakage errors 25], i.e., unwanted transitions 
to states that are outside the subspace spanned by the 
logical qubit states |0) and |1). The description presented 
here is an extension of earlier results on the SC flux qubits 
p2i and has potential applications to hybrid charge-flux 
qubits 7]. The role of the self and mutual inductances 
in SC charge qubits have been previously studied [2^ , in 



particular as a means of coupling two SC charge qubits P, 
|24| . Here, we fully and systematically take into account 
self and mutual inductances in the underlying SC circuit. 

While the circuit theory developed in Secs. HlllHIl and 
1^3 can be applied to any SC charge qubit, its usefulness 
will be illustrated with some specific examples of charge 
qubit circuits that have been studied before in Sec. 
where we reproduce and extend some previously known 
results. However, we stress that the circuit theory results 
are more general than previously applied methods for the 
following reasons, (i) The derived Hamiltonian is not a 
priori truncated to a two-dimensional subspace, which al- 
lows us to treat leakage and to derive the matrix element 
of the system-bath coupling, (ii) The capacitance matrix 
of the circuit is fully taken into account, and no assump- 
tion about the relative magnitude of gate and Josephson 
capacitances has to be made, (iii) The inductance matrix 
of the circuit is fully taken into account. 

Any number of dissipative elements Z (external 
impedances, resistances) can be included in the circuit 
theory. In our treatment of the system-bath Hamilto- 
nian and the decoherence and relaxation rates in Sec. lIVI 
we choose to restrict ourselves to the case of a single 
impedance Z in order to keep the notation simple. How- 
ever, the analysis can readily be extended to multiple 
impedances in analogy to SC flux qubits |27| . 
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FIG. 1: Circuit graph of a single voltage-biased charge box. 
Branches represent a Josephson junction (Ej), capacitances 
(Cj and Cg), a voltage source V, and the impedance Z. The 
nodes are shown as black dots; the node connecting the junc- 
tion (Ej) to the gate capacitance Cg represents the SC island. 
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II. NETWORK GRAPH THEORY 

The purpose of this Section is to derive KirchhofT's laws 
and the current-voltage relations (CVRs) for the circuit 
of a general SC charge qubit in an appropriate form for 
their later use in the derivation of the classical equations 
of motion of the circuit fSec. 

Our analysis (see also |2^]) starts with the represen- 
tation of the SC circuit as a directed graph, in which 
the branches 61 , 62 , • • • , represent one of the follow- 
ing lumped circuit elements: a Josephson junction, in- 
ductance, capacitance, voltage source, or an external 
impedance (e.g., a resistance). The circuit graph of a 
single, voltage-biased charge-box in Fig.^is a simple ex- 
ample of a circuit graph. In our examples, we neglect the 
quasiparticle resistance of the junctions because it makes 
the analysis simpler and because they are typically less 
important than the impedances of the voltage sources; 
however, the shunt resistances can easily be included as 
additional impedances in the circuit. The next step is 
to find a tree of the graph, i.e., a loop- free subgraph 
which connects all nodes (for each connected piece of the 
graph, we choose a connected sub-tree). The branches 
/i> /2j • ■ • , /f outside the tree are the so-called chords; 
each chord /i, when added to the tree, gives rise to a 
unique loop, a fundamental loop Ti of the circuit. The 
topological information about the graph which is of im- 
portance for our analysis can be represented in the fun- 
damental loop matrix (i = 1, . . . _F; = 1, . . . , _B), 



1, if hj E Ti (same direction), 
— 1, \ihj E Ti (opposite direction), 
0, 



(1) 



where the direction of the fundamental loop Ti is defined 
to be opposite to the direction of its defining chord /i. 
Accordingly, the currents I = {Ii, . . . , Is) and voltages 
V = (Vi, . . . ,Vb) associated with the branches I, . . . ,B 
of the graph are split up into tree and chord currents and 
voltages. 



I=(Itr,Ich), 



V= (Vtr,Veh). 



(2) 



With this ordering, the fundamental loop matrix assumes 
the form 



F(^) = (-F^ 1 1) , 



(3) 



and we will simply refer to the matrix F in the following. 
Using Eq. ij^Lj^s write Kirchhoff's laws in the following 
useful form 



FIch — ^Itr, 
F^Vtr = Vch-*., 



(4) 
(5) 



where = (*&!: ■ • ■ j ^f) denote the externally ap- 
plied magnetic fluxes threading loops 1, . . . , F of the cir- 
cuit. The partition of branch types into tree and chord 
branches is dual to the flux qubit case i.e., the roles 
of tree and chord branches are interchanged. 



Before we proceed, we summarize the assumptions 
about the circuit that will be used in the following. 

(i) There are no loops containing nothing else than 
Josephson junctions (J), external impedances (Z), 
and voltage sources (V). This assumption is physi- 
cally motivated because all loops have a finite self- 
inductance. 

(ii) Voltage sources (V) and impedances (Z) are not 
inductively shunted. 

(iii) There are sufRciently many capacitors (C) in the 
circuit to independently shunt all inductors. A 
more precise form of this requirement is that the 
capacitance matrix C has full rank (see below). 

Using assumption (i), we may split up the current and 
voltage vectors as 

It,. = (Ij,lL,Iy,Iz), Ich = (Ic.,Ic,I/f), (6) 

Vt,. = (Vj,VL,Vy,Vz), Veh = (Vc,,Vc,VK). (7) 

The chord current and voltage vectors Ich and Vch in 
Eqs. jnj and Q contain the branch currents and volt- 
ages of the capacitors (Cj,C) and chord inductors (K); 
the tree current and voltage vectors Itr and Vtr contain 
the branch currents and voltages of tree inductors (L), 
Josephson junctions (J), external impedances (Z), and 
bias voltage sources (v) jl^. The loop matrix F then 
acquires the block form. 



/ 1 Fjc FjK 

Flc Flk 

Fyc FvK 

\ Fzc FzK 



(8) 



The form of the first column in Eq. (jS]) reflects the fact 
that the Cj capacitances are (by definition) shunted in 
parallel to the Josephson junctions. Moreover, assump- 
tion (ii) above imphes Fvk — Fzk = 0. In order to 
derive the equations of motion, we formally define the 
branch charges and fluxes {X — C,K,J,L, Z, V), 



lx{t) - Qx{t), 



(9) 
(10) 



where the formal fluxes of the Josephson branches are 
the SC phase differences across the junctions, according 
to the second Josephson relation. 



^.7 ^ y 
$0 27r' 



(11) 



with <I>o = h/2e. The current-voltage relations (CVRs) 
of the Josephson, capacitance, and external impedance 
branches are 



Ic sinKf = If. sin 2tt 



Qc = CVc, 
Vz = Z*Iz, 



(12) 

(13) 
(14) 
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where the convolution is defined as (f *g)(t) = /Iqo f(i~ 
T)g{T)dT. The CVR for the inductive branches has the 
following matrix form, 



II 
Ik 



1/ 



, (15) 



where L and Jjk are the self inductances of the chord 
and tree branch inductors, resp., off-diagonal elements 
describing the mutual inductances among chord induc- 
tors and tree inductors separately, and Lilk is the mu- 
tual inductance matrix between tree and chord inductors. 
Since the total inductance matrix is symmetric and pos- 
itive, i.e. v-'^Ltv > for all real vectors v, its inverse 
exists, and we find 



Ik 



t-1tT f-1 
'^K ^LK^ 



■'K 



with the definitions 



L — L — Li/f L^^L"^ 



^K — 



La L 



LK 



LK^ 

L-iL 



LK- 



(16) 



(17) 
(18) 



III. CLASSICAL EQUATION OF MOTION 

In this Section, we derive the classical equation of mo- 
tion of the dynamical variables $ = ($,/,^'i) of the 
circuit. 

We now combine Kirchhoff's laws, Eqs. I0J and lO, 
and the the CVRs, Eqs. lfT^ - ((TH|l . in order to derive the 
classical equations of motion of the circuit. These will 
then be used in Sec. IIVI to find the Hamiltonian suit- 
able for quantization. The details of the derivation are 
explained in Appendix fXl 

Equations (|A3|) and (jA9p can be summarized as 

C* = Q-CyV- J^cCz*Vc, (19) 

with the combined flux vector # = ($,7,^>i) ~ 
(<f>o¥'/27r, *l), and the canonical charge 



Q = - 



Qj 

Ql 



K- 



(20) 



Note that in the SC charge qubits studied in Ref. ^ the 
Josephson junctions lead to (otherwise only capacitively 
coupled) SC islands, with the consequence that there are 
no chord inductors (K), and Q — — (Qj, Ql)^. However, 
the quantronium circuits which have hybrid charge 
and flux nature, cannot be described without chord in- 
ductors. In the following, we will derive our theory for 
the most general case including chord inductors, but fur- 
ther below, we will also discuss the much simpler special 



case without chord inductors. In Eqs. H19|l and H20() we 
have also introduced the notation 



JX 



for X — C, K , and the capacitance matrices 
Ctot CjL \ _ I u \ , ^ 



(21) 



C 



Ij+TcCT^, (22) 

(23) 
(24) 



We can further rewrite the dissipative term in Eq. H19() 
by using Eq. © (capacitance part), solving for Vp, and 
substituting the solution back into Eq. (|19|l . with the 
result 



{C + Cz)*1> = Q-CvV, 
where we have introduced 



Cziio) = mCz(t^)m^, 



rh = TcCF 



zc 



m,7 



(25) 

(26) 
(27) 



Cz{uj) = iujZ{Lu) (l + FzcCFlcii^Ziuj)) .(28) 

Using the symmetry of Cz^i^), we can show that Cz{i^) 
is also a symmetric matrix. 

We obtain the equation of motion from Eq. 125|) by tak- 
ing the derivative with respect to time, and using Eq. Q 
with X = K,L and Eq. ((TS|l . 



dU 

(C + C.)** = Q = --, 



(29) 



with the potential 

[/(*) = -L7^cosv3 + i*^Mo* + *'^N*^, (30) 

where $ = (<i>oy/27r, ^•l) and 

Mo^GL^^g^, N = GL7'{0 1a )^, (31) 

with the {Nl + Nk) x {Nj + Nl) block matrix 

-FjK 



G = 



1l -Flk 



(32) 



Using L( = Li, we observe that Mq = Mo. In the 
absence of chord inductors (K), we find ^-^Mq^" = 
#jL~^#i and N = 0, whereas in the absence of 
tree inductors (L), we obtain i$^Mo$ + ^^N^^ — 

1 (FJa-'P + *.)'^ La (f5k<P + *.) + const. 

By bringing the dissipative term in Eq. ^ to the 
right hand side and using assumption (iii), we find the 
equation of motion 



dU 



(33) 
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with the dissipation matrix 

= (1 + C2HC-i)"'CzH = mKHm^, (34) 
and the frequency-dependent kernel 

K{uj) = Cz{(^) (l + m^C-^mCz{uj)y^ . (35) 

Since both 0^(0;) and C are symmetric matrices, we find 
that K(a;), and thus also Cdii^), are symmetric. More- 
over, we know that Cd(t) inherits two additional proper- 
ties from Z{t): it is also real and causal, i.e., Cd{t) = 
for t < 0. In a perturbation expansion in Z^, the lowest 
order term in K(ti;) is simply K((jj) = iujZ{u!) + 0(Z)^. 

In deriving Eq. H33(l , we have used assumption (iii) that 
the matrix C has full rank, such that exists. Since all 
junctions are capacitively shunted, we know that Ctot has 
full rank, hence Nj < rankC < Nj+Nl, where Nx is the 
number of branches of type X. The case rankC < Nj + 
Nl occurs if there are not sufficiently many capacitances 
in the circuit to independently shunt all inductors. In 
that case, Eq. (|33|) . without the dissipative part, contains 
I = Nj + Nl — rankC constraints that can be used to 
eliminate / degrees of freedom. In the case of SC fiux 
qubits [i^l, it was assumed that only the junctions are 
shunted by capacitors (rankC = iVj), thus I is the number 
of tree inductors. 

IV. QUANTUM THEORY 

The purpose of this section is to derive the Hamil- 
tonian of the circuit, including its dissipative elements, 
and then to quantize this Hamiltonian in order to have 
a description of the quantum dissipative dynamics of the 
circuit from which a master equation and, finally, the 
decoherence rates can be derived. 

The Hamiltonian of the circuit 



1 



- ^ (Q - CvYf (Q - Cy V) + [/(*), (36) 

giving rise to the equation of motion (|33|l without dissi- 
pation (Z = 0), can readily be quantized with the com- 
mutator rule 



Qj] = 



I] ■ 



(37) 



A somewhat subtle point here is that while the inductor 
flux variables #i are defined on an infinite domain, the 
Josephson flux variables #j = (<I>o/27r)y> are deflned on 
a compact domain since they are periodic with period 
^o- Upon imposing Eq. (|37|l . this leads to charge op- 
erators Ql with a continuous spectrum and Qj with a 
discrete spectrum with eigenvalues Qji — 2eni, with rii 
integer 

In order to describe the dissipative dynamics of the 
SC circuit, we construct a Caldeira-Leggett Hamiltonian 
po] H = Hs +1-13 + "Hs B that reproduces the classical 
dissipative equation of motion Eq. (|33|) . For simplicity. 



we will restrict ourselves to the case of a single impedance 
Z here, where a single bath of harmonic oscillators can 
be used to model the dissipative environment. 



Pa 



1 



2^^ + 2™""'^' 



(38) 



We choose the system-bath coupling to be of the form 

T-i-SB — C^m ■ Q ^ CaXa = m ■ C^^Q ^ CaXa, (39) 
a a 

such that it reproduces the classical equation of motion 
Eq. (|33l) . with a spectral density of the bath modes (for 
a derivation, see Appendix Ib)) 



J{uj) = -1tjiK{uj). 



(40) 



Note that the kernel K has become a scalar because we 
are now only dealing with a single external impedance. 

From the Hamiltonian 7i, the master equation for the 
evolution of the system density matrix can be derived 
|22| . In the Born-Markov approximation, the matrix el- 
ements pnm = (n\pf;\rn), where Hsln) — ujn\n), obey the 
Redfleld equation jSfl 



{t)~Y,^^rnklPkl{t), (41) 



kl 



with I 



~ LOn — t^m, and with the Redfleld tensor, 



Rnmkl — ^Ira ^ ^ nrrk ^ ^nk ^ ^ F 



.(-) 
Irrm 



.(+) 

Imnk 



F 



where {r[ti,r 



(-) 
Imnk ' 

(42) 



ReF 



(+) _ 



-(m • Q)/„i(m • Q)„fe J(|w„fc| 
n 



ImF 



(+) _ 



Imnk 



= --(m • Q);„j(m • Q)„fe X 



P du- 



UJ — OJnk COth 



sinh?i/3|w„fc|/2 ' 
(43) 



and m = C^^m. 

The Redfleld equation l|41() can be derived for arbitrary 
SC circuits. The SC circuit can represent a single qubit 
or a number of qubits. In order to make connection with 
single-qubit experiments, we apply the theory to the case 
of a SC circuit representing a single qubit. Restricting 
ourselves to the two lowest levels and working in the sec- 
ular approximation the Redfleld equation Eq. |(1T)) 
turns into a Bloch equation with the relaxation (Ti) and 
decoherence (T2) times. 



1 

1 _ 1 



|(0|m-Q|l)|V(woi)coth 



1 

t1 



^ = i|(0|m.Q|0)-(l|m.Q|l)p 

Id) Ti Tllu 



2kBT' 



J{co) 



(44) 
(45) 



2kBT. (46) 
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In the semiclassical approximation |23|, (0|Q|1) « 
(l/2)(A/c^oi)AQ and (0|Q|0) - (1|Q|1) « (eMi)AQ, 
where AQ = Qo — Qi is the "distance" between two 
locahzed low-energy classical charge states Qo and Qi, 
e is the classical energy difference and A the tunneling 
amplitude between them, and ujqi = VA^ + is the en- 
ergy splitting between the two quantum eigenstates in 
this energy double well. Within this approximation, we 
find 



1 
1 



m AQ|2 
m AQ|2 



A 

e 



J(ti;oi) coth 



TlLUoi 



2fcsT 



huj 



2fcsT. 



(47) 



(48) 



The leakage rates from the logical state fc = 0, 1 to states 
n = 2, 3, . . . outside the computational subspace can be 
estimated as 



f 



-J2\(k\m-Q\n)\\jMcoth^. 

n=2 ^ 



(49) 



V. EXAMPLES 

A. Single Charge Box 

The voltage-biased charge box is shown in Fig.^ where 
the inductance of the leads has been neglected for sim- 
plicity (no L and K branches). The tree of the graph is 
given by the Josephson, voltage source, and impedance 
branches. For the loop matrices, we simply find 



Fjc = Fyc = F 
With the capacitances 

C = Ctot = Cj + Cg, 

we arrive at the Hamiltonian, 

{Qj+CgVf 



zc 



Us = 



2a, 



Cv — Co 



Ej cos (p. 



(50) 



(51) 



(52) 



The coupling to the environment is characterized by m = 
(Cg/Ctot)- As an example, we give here the relaxation 
and dephasing times, with m — |m| — Cg/Ctot, 

— = 2nm 4|(0|n|l)| t^o^ coth — — , (53) 



1 

t1 



2W|(0|n|0) 



where n = Q/2e and Rq — h/e^. In the semiclassical 
limit, (0|n|l) (l/2)(A/woi)An and (0|n|0) - (l|n|l) « 
(e/a;oi)An. With An k, 1, we reproduce the results in Pi- 
Typical leakage rates are of the form of l/Ti, with the 
matrix element replaced by |(0|n|fc)| and |(l|n|fc)|, where 
k > 2 labels a state other than the two qubit states, and 
with loqi replaced by ujik {I = 0, 1). 




FIG. 2: A flux-controlled Josephson junction. 





V, 



I I 



X 



FIG. 3: Two inductively coupled charge boxes. 



B. Flux-controlled Josephson junction 

A flux-controlled Josephson junction is a SC loop with 
two junctions which acts as an effective Josephson junc- 
tion with a flux-dependent Josephson energy |5l|. The 
circuit Fig. [3 we use to describe the the flux-controlled 
junction comprises a chord inductance (K) with induc- 
tance L. The tree consists of the two Josephson branches. 
The only relevant loop matrix is FjK = ( 1 -1 ) .In 
the hmit L 0, and if Eji = Ej2, we find F'^j^ip + ^x = 
ifi ~ ip2 + ^ 0, which leads us to the Hamiltonian 



■Hs 



91 

2C 



E J {<^x) cos ip, 



(55) 



where (p — ipi + 7r$j;/$0j 
2Ejcos{2tt<Px/'^o)- 



C^C.n+Cj2, and E.j{^x) 



C. Inductively coupled charge boxes 

We now turn to the case of two charge boxes of the 
type discussed in Sec. IV Al coupled via an inductive loop 

nil, 

as shown in Fig. |2| Here, the tree consists of all 
Josephson, voltage source, and impedance branches, plus 
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the inductive branch L, and the loop matrices are 
Fjc = Fyc = Fzc = J 5 ) , Flc = ( 1 1 ) . (56) 

With the two capacitance matrices C = diag(Ci, C2) and 
Cj = diag(Cji,Cj2), we find Ctot - C + Cj, Cjy = C, 

I 



CjL = Cly = (Ci,C2)^, and = Ci + C2. The 
vector m consists of the two parts mj = C and 111^ = 
( Ci C2 ). With Eq. ^ and the inverse of the total 
capacitance matrix, 



1 

7 



C2)Cj2 
C1C2 
— CiCj2 



CIC2 

C2)Cji 

— C2CJ1 



— Ci Cj2 
CjlCj2 



'-'cff,! 
'-'cff. 12 



'-'cfT,12 
"-"off, 2 



'-'cfT.Ll 
'-'cfT,L2 



(57) 



'-'cff, LI '-^cff,L2 '-'cff.L 



where 7 = {Ci+C2)CjiCj2-ClCj2-ClCju the Hamil- 
tonian of the coupled system can be written as, 

= 2^ 1 — \- Ejt cos (p, 



\ 2Ccff i 

{Ql + CiV, + C2V2? ^ 

2CcfT,L 2L 

(Qjl+Cl'^l)(Qj2+g2^2) 



(58) 



Ccffa2 

iQ.n + a^^)(QL + CiVi + ^2^2) 

Ccff,Li 



■E 

i=1.2 



While the last term in Eq. H58|l couples each qubit to 
the LC mode associated with the inductor L, and is thus 
responsible for the inductive coupling of the qubits, the 
second last term provides a direct capacitive coupling 
between the qubits. In the limit Ci <C Cji, we reproduce 
the results of 1]; however, there are additional terms 
of order Ci/Gji, in particular the new term (x l/Coft,i2 
in the Hamiltonian that capacitively couples the qubits 
directly. Since the coupled system involves at least four 
levels (more if excited states of the LC coupling circuit 
or higher qubit levels are included), it can no longer be 
described by a two-level Bloch equation with parameters 
Ti and T2. We can however fix one of the qubits to 
be in a particular state, say |0), and then look at the 
"decoherence rates" of the other qubit. To lowest order 
in Ci/Cji, these rates due to the impedance Zi have the 
form (g, = C,/(Ci + C2)) 



1 
1 



APPENDIX A: DERIVATION OF THE 
EQUATIONS OF MOTION 

This appendix contains the derivation of Eq. H19|l . 
Note, first, that the externally applied magnetic flux 
only threads loops with a finite self- inductance (i.e., those 
pertaining to a chord inductor, K), and not, e.g., the cir- 
cuit loop formed by a junction J and its junction capac- 
itance Cj, therefore x = (*cj'*C'*if) = (OiO,$|-). 
Using this fact and Eqs. jsj (capacitance part) and 
we obtain 

I^FJc^ = Vc - F^cVl - F^pVy - ^Ic^z (Al) 
zvr 



Rq 
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C iQ 



c - Flc*L - F^cVy - F'^cZ * Iz, 



multiply this equation by FjcC and use Eq. Q 
(impedance part), with the result 



$0 
27r 



-FjcCFyj^Vy 



CFicZFzc * Qc- (A2) 
Then we make use of Eq. Q (Josephson part) and obtain 

$0, 



2tt 



— C./vVv — FjcCz * Vf7, 



(A3) 



where we have defined the frequency-dependent capacity 
Cz{uj) = iujCFl^Z{Lu)FzcC and 



2vrgf4|(00|n,|10)p '^^y""^) c.oicoth|^, (59) 



(60) 



Ctot = Cj + FjcCFJ(^, 

CjL = FjcCF 
C,7v — FjcCF 



vc^ 



(A4) 
(A5) 
(A6) 



We find that Cz{'jj) is a symmetric matrix since both C 
and Z are symmetric. Using Eq. ^ (capacitance part) 
again, we obtain 



Flc^L = C-iQc - 



2tt 



JC'P 



'vc^v 



(A7) 
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which we multiply with FicC, with the result 

FlcCFIc^l = FlcQc-^FlcCF^c'P 

Ztt 

—FlcCFyuVv 
-FicCFl^ZFzc * Qc- (A8) 

With the definitions Cl = FlcCF]]^ and Clv = 
FicCF^c we find 



Cl*l = -Ql-FlkQk -^C^^if 



(A9) 



Equations IjASp and (|A9|I are rewritten in a more compact 
form in Eq. ((T^ . 



APPENDIX B: SYSTEM-BATH DYNAMICS 

In this section, the form of the system-bath coupling 
operator Hsb and its spectral density J{lu), Eqs. 
and H4U|I . are derived in detail. 

We first inspect the Hamilton equations for the bath 
coordinates, 



OH _ Pa 



(Bl) 



Pa = = -maUjlxa - CaVh- C^^Q, (B2) 

then take their derivative with respect to time, and solve 
them in Fourier space. We obtain 



Xaiuj) 



CqIii-C 

ma [lo'^ - Wq) ' 



(B3) 



iujCam ■ C 

Pa[UJj = malLUXa[UJ) = 5 5 ■ (B4) 



Next, we look at the Hamilton equations for the system 
coordinates, 

* = 1^ =C-i |^Q + m^c„a;„j , (B5) 

(B6) 



dH 


dU 


"a* 




. 113 


and llB6l 



we obtain 



C* - + - - if * m(m • C-iQ) 



nia dip 



where 



(B7) 



(B8) 



directly determines the bath spectral density 



J(Lo)^^y"-^S(uj-uJa) = -lmK(u;). (B9) 

2 '^-^ TTlaUJa 
a 

By comparing Eq. HB7|) with Eq. (123, we find 

Cd{uj) = if (w)mm^. (BIO) 
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